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Abstract. Using the observation that configurations of iV polymers with hard 
core interactions on a closed random surface correspond to random surfaces with N 
boundary components we calculate the free energy of a gas of polymers interacting 
with fully quantized two-dimensional gravity. We derive the equation of state for 
the polymer gas and find that all the virial coefficients beyond the second one vanish 
identically 
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It is well known that some statistical mechanical models are easier to solve on 
a random triangulated surface than on a fixed lattice. Well known examples are 
the Ising model |lj and self avoiding random walks ||. In the latter reference 
a closed form for the partition function of n mutually self avoiding walks with 
common endpoints (watermelon diagrams) was obtained and the critical exponents 
determined; see [pj for a review. 

The reason for the simplicity of self avoiding random walks on a random surface 
is that we can keep the walks fixed and sum over the surfaces in a relatively straight- 
forward fashion. One can for example regard the links which make up the walks as 
boundary components where the links have been identified pairwise. In fact one can 
also identify the boundary links in such a way that one obtains a branched polymer 
on the surface. The statistical sum over surfaces and walks therefore becomes a sum 
over surfaces with fluctuating boundaries. 

In [U closed expressions were derived for all the multiloop functions for random 
surfaces made up of even sided polygons. The scaling limit of the loop functions 
in these models was also calculated and it was shown that the limit is the same 
when one drops the assumption that the polygons are even sided. The authors 
of [|j] used matrix model techniques. Identical results can be obtained working 
directly with triangulations ||. Here we use the explicit results of H to study 
the thermodynamical properties of a gas of polymers interacting with 2d quantum 
gravity. This might be regarded as a two-dimensional toy model for primordial 
cosmic strings in thermal equilibrium with the gravitational field. 

Let us consider an ensemble of triangulated surfaces with the topology of S 2 with 
b boundary components. By a triangulation we here mean a collection of polygons 
whose sides are identified pairwise and degenerate boundaries may appear (unre- 
stricted triangulations in the terminology of 0). We assume that each boundary 
component consists of an even number of links and carries one marked link. Let 
us arrange the links cyclically along each boundary component: ii, i 2 , . . .i^i-, so that 
%\ is the marked link. If we identify the link %\ with 121, 12 with 121-1 etc. we ob- 
tain a self- avoiding walk on a surface with one less boundary component. One of 
the enpoints of the walk is marked corresponding to the marking of the original 
boundary component. Doing this for all the boundary components we end up with 
a collection of b mutually self-avoiding walks. We shall call the self-avoiding walks 
linear polymers or just polymers from now on. Given an orientation of the surface 
we see that there are precisely two distinct markings of a boundary component that 
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give rise to the same unmarked linear polymer. Similarly one can obtain a branched 
polymer on the surface by identifying the boundary links in a more general fashion. 
This can be done in precisely 

r W m 
Cl ~ WTTy. (1) 

different ways, see e.g. ||. Here we shall focus on the case of linear polymers. The 
general case can be treated by the same methods and we give the result below. 

It follows from the above considerations that we can write the partition function 
of a gas consisting of b polymers of lengths h,h,...,lb with a hard core interaction 
on a random surface as 

Z b (h,...J b ) = 2- b w(2h,...,2l b ), (2) 

where w is the 6-loop function of quantum gravity (discretized Hartle-Hawking wave 
functional). The loop function is given by 

w(l u l 2 ,...,l b ) = Yl e ~ ST ( 3 ) 
TeT(h,...i b ) 

where T(h, . . . , l b ) is the collection of all triangulations in the class under study 
which have b boundary components of lengths l±, . . . , l b and St is the action of the 
triangulation T. 

We briefly recall the derivation of the formula for the loop function. For simplicity 
we consider triangulations which are made up of even sided polygons and take the 
action to be of the form 

oo 

S T = -Y l k j ]ogg 2j , (4) 

i=i 

where kj is the number of 2j-gons in T. Only a finite number of the coupling 
constants g 2 j are nonzero and they are assumed to be nonnegative since we are 
interested in the pure gravity case but in principle one can study the multicritical 
models by the same method. Denoting the coupling constants collectively by g we 
can write the Laplace-transform of the 6-loop function as 

b d 

w{g\ z u . . . , z h ) = J] (5) 

i=2 aV \ Zi > 

where 

oo 

<%*) = $>- 2 <-Vs;20 (6) 
i=i 
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is the Laplace-transform of the one-loop function and the "loop insertion operator" 
is defined by 

d 00 d 

see j|, [|. One should regard the variables Z{ as the cosmological constants for 
boundaries. After some calculation one obtains for b > 3 the formula 

a{9 - :zi -' Zb) = (m^i) s^n^ (8) 

where c = c(g) is the critical value of the Zj- variables, 

2 / duj ujV'{uj) 



and 

00 

V{u) = Y.92i" 2j . (10) 

i=i 

Writing g 2 j = gt 2 j where g > and the t 2 j are kept fixed we have a critical value go 
of g where the loop functions become singular. 

The scaling limit is calculated by introducing a lattice length a which tends to 
zero and letting g ^ go and c — > Co = c(go) such that 

log#-log# = o 2 A (11) 
log Zi- log c = aZ h (12) 

where A and Zi are renormalized dimensionfull cosmological constants which are 
kept fixed in the scaling limit. The continuum loop functions are given by 

W(A;Z u ...,Z b ) = \ima- 5+n / 2 w(g;z 1 ,...,z b ) 



d \ b ~ 3 1 A 1 



" C [dK) VxUi{z i + VX)W' (13) 

for b > 3 and C is a numerical constant. There are also explicit formulas for 6=1 
and 2 but they are not of interest to us here. 

It is now straightforward to take the inverse Laplace-transform in the variables 
Zi and A and this yields the continuum 6-loop function at fixed volume V and loop 
lengths L\, . . . , L b : 

W(V; Li,...,L b ) = CV h ' 7 / 2 f[ ^ exp (~ ^ L^j \ . (14) 
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Let now all the polymers have the same continuum length L. Then the partition 
function becomes 

T/6-7/2 

Z b (V,L) = C^—— L"/ 2 ^ 2 /^, (15) 
bi 

where we have inserted the statistics factor in order to ensure "correct Boltz- 

mann counting" of polymers. Putting 

P=\ (16) 

and taking the limit b — > oo and V — > oo with the polymer density p fixed we obtain 
the free energy per unit volume 

/(p,L) = plogp-^logL + ^. (17) 

The first term on the right hand side of (|H]) corresponds to an ideal gas. 
The pressure is given by 

d 



n 2 T 2 

= P+ Pj ^- (18) 

The equation ([18]) is the exact equation of state for the polymer gas and we see in 
particular that the hard core repulsion of polymers is reflected in a single quadratic 
term in the virial expansion of the pressure. This is much simpler than the cor- 
responding formula for the pressure for polymers in flat space where in general all 
the virial coefficients are non- vanishing |6|, [7]] . It is clear that this gas has no phase 
transition as expected for a gas with purely repulsive interaction. 

In the case of general branched polymers we obtain an extra entropy factor 



(19) 

on the right hand side of (0). Using 

d ~ 4'r 3/2 (20) 

for large I we find that the continuum partition function at fixed volume is in this 
case given by 

T/6-7/2 

Z' b (V, L) = C—y—L- b e- b2L2/4V , (21) 
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which yields the free energy 

2 t 2 

f'(p, L) = p log p + p log L + P — (22) 

and the equation of state is unchanged. 

Using the semi-explicit formulas for the higher genus loop functions derived in 
|| one can in principle generalize the above calculation to obtain the free energy 
of a polymer gas on a surface of arbitrary fixed genus. It is not straightforward 
to obtain an explicit expression in this case since there is no known analogue of 
(P|) for higher genus surfaces. The problem of letting the number of handles go to 
infinity in the thermodynamic limit so that one has a finite density of "wormholes" 
can presumably also be studied by this method. 
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